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ABSTRACT 

An increasing amount of observational evidence supports the notion that there are two modes of star formation: a quiescent mode 
in disk-like galaxies, and a starburst mode, which is generally interpreted as driven by merging. Using a semi-analytic model 
of galaxy formation, we derive the relative contribution to the cosmic star formation rate density of quiescently starforming and 
starburst galaxies, predicted under the assumption that starburst events are triggered by galaxy encounters (merging and fly-by 
kind) during their merging histories. We show that, within this framework, quiescently starforming galaxies dominate the cosmic 
star formation rate density at all redshifts. The contribution of the burst-dominated starforming galaxies increases with redshift, 
rising from < 5 % at low redshift (z <0.1) to ~ 20% at z > 5. We estimated that the fraction of the final (z=0) galaxy stellar mass 
which is formed through the burst component of star formation is ~ 10% for 10 10 M o <M» < 10"- 5 M o . Starburst galaxies, selected 
according to their distance from the galaxy main sequence, account for ~10% of the star formation rate density in the redshift 
interval 1.5< z <2.5, i.e. at the cosmic peak of the star formation activity. 

Key words, galaxies: evolution - galaxies: fundamental parameters - galaxies: interactions - galaxies: starburst 



1. Introduction 

How galaxies build up their stellar mass is a central ques- 
tion in galaxy formation. Many previous works have demon- 
strated that a substantial fraction (30%-50%) of the stellar mass 
formed between z ~ 3 and z ~ 1 (Dicki nson et al.||2003 



Fontana et al.|20031|2004l|2006l|Glazebrook et al 
et al 



2004; 



IRudnick et al.J|2006| 



Papovich et a 



2004l|Drory | 
. 2006TTYlffl1 



et al.|f2006[ [Pozzetti et al.||2Q07) . However, it is still unclear 



what is the major channel for galaxy growth during this pe 
riod. Two star formation modes have been identified: the first 
is a "quiescent" Qmode in disk-like galaxies where the star for- 
mation is extended over the whole galaxy disk and occurs on 
time scales of ~ (1 — 2) Gyr ; while the second is a compact, 
starburst mode, where the star formation tends to be dominated 
by the nuclear regions and the gas depletion time scales are 
significantly lower (~ 10 7 - 10 8 yr) than those in quiescently 
starforming galaxies. This accelerated mode of star formation 
is likely triggered by major mergers, as suggested by obser- 
vational evidence and theoretical arguments. Observationally, 
major mergers are associated with enhancements in star for- 
mation in local Ultra Luminous Infrared Galaxies (ULIRGs, 
|Sanders & Mirabel fl996[|Elbaz et al.|2007|l, and some submil- 
limeter galaxies (SMGs |Tacconi et al.|2008||Daddi et al.|2007 



|2009| |Engel et al.|2010) > . On the theoretical side, the merger- 
driven scenario for starburst galaxies is supported by the re- 
sults of high-resolution numerical simulations, which showed 
the effectiveness of galaxy major mergers in causing the loss of 
angular momentum of the gas in the galaxy disk, with the con- 
sequent trigger of starburst events into the central regions of 
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1 In this paper, we use the term quiescent to refer to starforming 
galaxies that are not in interaction. We do not mean non-starforming 
systems, as the term is used in some literature. 



the galaxy (|Hernquist|[T 989 ; Bar nes & Hernquist|199i1 [1996] 



Mihos & Hernquist 1994 1996); and by the success of semi- 
analytic models (SAMs) of galaxy formation, which include 
this additional channel of star formation, in reproducing the 
early formation of stars in massive galaxies (e.g. Somerville 
|et al.|2001l|Menci et al.|2004||20fJ5l l. 

An effective tools to understand the relative contribution of the 
different star formation modes is provided by the scaling re- 
lations connecting the star formation rate (SFR) with global 
galactic quantities, like the gas mass (the Schmidt-Kennicutt 
relation) or the total stellar mass (M,). In fact, the normal- 
ization and/or the scaling of the SFR with such quantities dif- 
fers for starbursts and quiescently starforming galaxies. As for 
the former, the extension of the well-established local relation 
( |Kennicutt|T 998 ) to higher redshifts has shown that quiescently 
starforming galaxies follow the local Schmidt-Kennicutt rela- 
tion, while the SFRs of starburst galaxies are typically one or- 
der of magnitude greater than expected from their projected gas 
surface density ( |Daddi et al.|2010||Genzel et al.|2010) . As for 
the correlation between the SFR and M„ it has recently been 
shown that quiescently starforming galaxies lie along a "main 
sequence" characterized by a typical, redshift-dependent value 
of the specific star formation rate (SSFR=SFR/M„), while the 
less numerous starburst population is characterized by larger 
values of the SSFR. The main sequence is observed over a large 
range of redshift, from z ~ to z ~ 6, and it is characterized 
by a small scatter (~ 0.3 dex), by a normalization which in- 
creases by a factor of ~20 from z - to z — 2, and by a slope 
which is sensitive to the technique used to select the sample 
of starforming galaxies and to the procedure for measuring M* 
and SFR, with values ranging from 0.6 to 1 (Brinchmann et al. 
20041 |Noeske~et"aT. 2007; Elbaz et al |2007[|Daddi et al. 



2009; Santini et al. 2009; |Salim et al.|20071|Stark et al. 
|Gonz;ilezetal.|2011| >. " 
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The tecniques used to measure the SFRs and the stellar masses 
in galaxies are based on different observational tracers. A fruit- 
ful approach to estimate the stellar mass is the SED fitting. This 
technique is based on the comparison between a grid of spectral 
templates computed from standard spectral synthesis models 
and the available multiwavelength photometry (Papovich et al. 
2001] |Shapley et al.|[200T| [2005] |Fontana et al.||2006[ |Santini| 
et al.|2009] >. However, such a method does not provide reliable 
star formation histories at high redshifts, where the uncertain- 
ties become larger due to the SFR-age-metallicity degenera- 
cies. An alternative commonly used estimator of the SFR is 
the UV rest-frame band, where young and massive stars emit 
most of their light. Here the drawback is that dust absorbs, re- 
processes, and re-radiates UV photons at near-to-far IR wave- 
lengths, so that reliability of UV luminosity as a SFR tracer de- 
pends on large and uncertain corrections relying upon the dust 
properties. Since most of the energy radiated by newly formed 
stars is reprocessed by dust, in order to accurately derive the 
SFR it is necessary to determine the dust bolometric output. In 
the past, the total IR luminosity (from 8 to 1000 //m) was cal- 
culated by extrapolating observation in the mid-IR. With the 



launch of the Herschel Space Observatory (Pilbratt et al. 2010) 
it has now become possible to measure the far-IR luminosity 
of distant galaxies at wavelengths where the dust emission is 
known to peak. 

A recent study based on Herschel observations of starforming 
galaxies at 1.5< z <2.5 indicates that only ~2% of massive 
(M* > 10 10 M G ) galaxies in this sample have starburst nature 
and account for only ~10% of the star formation rate density 
(SFRD) at z ~ 2 ( |Rodighiero eTaT1[20TT] ). This finding rep- 
resent a test case for the merger-driven scenario for starburst 
galaxies, since in this framework the occurrence of the starburst 
events is driven by the galaxy merger rates which are expected 
to increase with redshift. In this paper we investigate this is- 
sue in the framework of a theoretical model of galaxy forma- 
tion that includes a physical description of starburst triggered 
by galaxy interactions during their merging histories. The lat- 
ter are described through Monte Carlo realizations, and are 
connected to gas processes and star formation using a semi- 
analytic model of galaxy formation in a cosmological frame- 
work HMenci et al.|2004l [2005| [20061 12008) . In the model the 
starburst events are triggered not only by galaxy major merg- 
ers (i.e. the fusion of galaxies with comparable mass, m * m') 
but also by minor mergers (m <K m') and by closer galaxy in- 
teractions which do not lead to bound merging (fly-by events). 
Minor mergers and fly-by events induce starburst with a lower 
efficiency compared to major mergers, but they are more proba- 
ble, so they contribute appreciably to the cosmic star formation 
history. 

A description of the SAM is given in Section [2] ; in Section [3] 
we derive the predicted SFR-M„ relation at 1.5< z <2.5 and 
the contribution to the cosmic SFRD of starbursts and quies- 
cently starforming galaxies; discussions and conclusions fol- 
low in Sections |4]and|5] 

2. The Model 

We u s e the SAM d escribed in details in |Menci et al.| ( |2004] 



2005 2006 



|2008[ ), which connects, within a cosmological 
framework, the baryonic processes (gas cooling, star forma- 
tion, supernova feedback) to the merging histories of the dark 
matter (DM) haloes. AGN activities, triggered by galaxy inter- 
actions in common DM haloes, and the related feedback pro- 
cesses are also included. Here we recall the basic points. 



2. 1. The dark matter merging trees 

Galaxy formation and evolution is driven by the collapse and 
growth of DM haloes, which originate from the gravitational 
instability of overdense regions in the primordial DM density 
field. This is taken to be a random, Gaussian density field with 
Cold Dark Matter (CDM) power spectrum within the "concor- 
dance cosmology" (Spergel et al. 2007) for which we adopt 
round parameters £2a = 0.7, Qq = 0.3, baryonic density 
= 0.04 and Hubble constant (in units of 100 km/s/Mpc) 
h = 0.7. The normalization of the spectrum is taken to be 
o"g = 0.9 in terms of the variance of the field smoothed over 
regions of 8 h Mpc. 

As cosmic time increases, larger and larger regions of the 
density field collapse, and eventually lead to the formation of 
groups and clusters of galaxies; previously formed, galactic 
size condensations are enclosed. The merging rates of the DM 
haloes are provided by the Extended Press & Schechter formal- 



ism (see |Bond et al.|[T991 Lacey & Cole 1993). The clumps 
included into larger DM haloes may survive as satellites, or 
merge to form larger galaxies due to binary aggregations, or 
coalesce into the central dominant galaxy due to dynamical 
friction; these processes take place over time scales that grow 
longer over cosmic time, so the number of satellite galaxies in- 
creases as the DM host haloes grow from groups to clusters 
(see |Menci et al.|2005||2006] l. 

2.2. The star formation law 

The processes connecting the baryonic components to the 
growing DM haloes (with mass m) are computed as follows. 
Starting from an initial amount mO^/Q of gas at virial temper- 
ature of the galactic haloes, we compute the mass m c of cold 
baryons within the cooling radius. The cooled gas mass m c set- 
tles into a rotationally supported disc with radius (typically 
ranging from 1 to 5 kpc), rotation velocity Vd and dynamical 
time td = rd/vd, all computed after |Mo et al. ( 1998] >. Two chan- 
nels of star formation may convert part of such gas into stars: 
i) quiescent star formation, corresponding to the gradual con- 
version of the cold gas in the galaxy disk into stars, for which 
we assume the canonical Schmidt-Kennicutt form: 



SFR q = m c /T Cj 



(1) 



where r q = qtd, and q is a model free parameter which is cho- 
sen to match the Kennicutt ( 1998 ) relation; 

ii) burst-like star formation triggered by galaxy interactions, 
at a rate: 



SFR b oc fm c /T b 



(2) 



here the time scale is assumed to be the crossing time for the 
destabilized cold gas component (fj), and / is the fraction of 
cold gas destabilized by the encounters. A clear visualization 
of the two star formation modes implemented in the model is 
given in figure [T] where we show the star formation histories 
(SFHs) for two subsets of galaxies selected from the Monte 
Carlo simulations. The fraction / of cold gas destabilized by 
the encounters has been worked out by Cavali ere & Vittorini| 
( 2000 ) in terms of the variation A j of the specific angular mo- 
mentum j as Gm/vd of the gas to read (Me nci et al.|20 04): 

I = 1 Vd \ 

2 I j I 2 \ m b V re i I 



f 



(3) 



here b is the impact parameter, evaluated as the greater of the 
radius r</ and the average distance of the galaxies in the halo, m' 
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Fig. 1. SFHs drawn from the Monte Carlo realizations of galaxies that, at z = , have masses of M« 1O 1O M (left) and 
M* 10 12 M Q (right). The SFHs are obtained by summing over all the progenitors that have merged to form galaxies of the 
selected final mass. The quiescent star formation is represented by the smooth component of the SFHs, the burst component by 
the impulsive one. The duration of the SFR-excess phase due to the interactions is represented by the width of the peak of the 
SFHs. 



is the mass of the partner galaxy in the interaction, V re i is the 
relative velocity between galaxies, and the average runs over 
the probability of finding such a galaxy in the same halo where 
the galaxy with mass m is located. The prefactor accounts for 
the probability 1/2 of inflow rather than outflow related to the 
sign of A j. We assume for the proportionality constant in eq. 
Q the value of 3/4, while the remaining fraction of the in- 
flow is assumed to feed the central black hole (see Sand ers &| 
Mirabel 1996). When applied to the black hole accretion and 
to the related AGN emission the above model has proven to be 
very successful in reproducing the observed properties of the 
AGN population from z = 6 to the present ( Menci et al.|2 004 
|2008l|Lamastra et al.|2010] i. 

The probability for a given galaxy to be in a burst phase is 
defined as the ratio T/,/T r of the duration of the burst to the av- 
erage time interval between bursts. The rate of such encounters 
TV 1 is 



r; 1 =n T I.( rt ,v,V rel )V rel (4) 

here rij = 3Nr/4nR 3 is the number density of galaxies in the 
same halo, V re i is their relative velocity, and Z ^ n(rj + r' t 2 ) is 
the cross section for such encounters which is given by Saslaw 
( 1985 ) in terms of the tidal radius r, associated to a galaxy with 
given circular velocity v (see Menci et al.|2004[ ). 



At each time step, the mass Am/, returned from the cold 
gas content of the disk to the hot gas phase due to the 
energy released by SNae following star formation is esti- 
mated from canonical energy balance arguments as Am/, = 
E SN €Qrj()Am t lvl, where Am, is the stellar mass formed in the 
time step, 770 is the number of SNe per unit solar mass (for a 
Salpeter initial mass function 770 = 6.5 x 10 _3 M Q _1 ), Esn = 
10 51 erg/s is the energy of ejecta of each SN, and v c is the cir- 
cular velocity of the galactic halo; eo = 0.01 is a tunable ef- 
ficiency for the coupling of the emitted energy with the cold 
interstellar medium. At each merging event, the masses of the 
different baryonic phases (Am c , Am* and Am/,) in each galaxy 
are refueled by those in the merging partners. Thus, for each 
galaxy the star formation is driven by the cooling rate of the hot 
gas and by the rate of refueling of the cold gas, which in turn is 
intimately related to the galaxy merging histories. The model 



has been tested against several observed properties of the lo- 
cal galaxy population like the B band luminosity function, the 
stellar mass function, the Tully-Fisher relation, the galaxy bi- 
modal colour distribution, the distribution of cold gas and disk 
sizes. At higher redshifts, the model has been tested against the 
evolution of the luminosity functions and the number counts in 
different bands ( |Menci et al.|2005l|2006] l. 



3. Results 

3. 1. The SFR-M* relation 

As a first step in the study of the contribution of starbursts 
and quiescently starforming galaxies to the cosmic SFRD, we 
compare in figure [2] the predicted SFR-M* (left) and SSFR- 
M, (right) relations for model galaxies with M, > 1O 9 5 M 
at 1.5 < z < 2.5 (coloured contours) with those obtained by 
Rodighiero etaL] pOTT) and by |Santini et"al1 ( p009] l. The for- 
mer authors used a combination of BzK- and Herscel-PACS- 



selected starforming galaxies from the COSMOS and GOODS 
fields. The BzK colour selection is designed to identify star- 
forming galaxies in the redshift interval 1.5 < z < 2.5 by using 
colours that samples key features in the spectral energy distri- 
butions of galaxies, mainly the rest-frame 4000 A break and the 
UV continu um slope (BzK= (z - K) AB -{B- z)ab >-0.2,[Daddi 
et al.|2004] l. The SFRs of BzK galaxies are estimated from the 



UV-rest frame luminosity corrected for dust reddening and they 
yield an almost linear correlation with the stellar mass (black 
circlces in fig. [2]). The Herschel-PACS detection limit depends 
on redshift; at z — 2 it corresponds to ~200 Ms/yr and to ~50 
M Q /yr for the COSMOS and GOODS fields respectively. The 
PACS-based SFRs run almost flat with the stellar mass and oc- 
cupy the upper envelop of the distribution (green squares in fig. 
|2j. This different behaviour could be ascribed to the different 
selection of the two samples, one being SFR limited, and one 
being mass limited (Rodighi ero et al.|201 l| l. The Santini et al 
(2009) sample is based on the GOODS-MUSIC catalog; for 
dectected sources the SFRs are derived from the 24/um emis- 
sion, corrected to take into account the well-known overesti- 
mation of the SFR for bright sources at z ~ 2 (see |Santini| 
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Fig. 2. Left: SFR-M, relation at 1.5 < z < 2.5. The five filled contours correspond to equally spaced values of the density 
(per Mpc 3 ) of model galaxies in logarithmic scale: from 10~ 6 for the lightest filled region to 10~ 2 for the darkest. Green squares: 
Herscel-PACS- selected galaxies, black circles: BzK galaxies, from the Rodighiero et al. ( 201 1 1 sample; blue triangles: GOODS- 



MUSIC catalog from the Sant ini et al.|(2009]) s ample. The dashed and long-dashed lines indicate the best-fit of the galaxy main 



sequence obtained by Rodighiero et al. ( 201 l| l and Santini et al. (2009) respectively, while the solid line shows the spline fit 



to the peaks of the SFR distributions of model galaxies as a function of stellar mass. Right: SSFR as a function of M». Colour 
code, symbols and solid line as in the left panel. The dashed and thin s olid lines indicate the SSFR > SSFRms-2ct and the 
SSFR> t H {zT l limits of the galaxy samples used in the analysis of Section 3.2 respectively. 



|et al.|2009] for further details), while for undetected objects the 
SFRs are derived from SED fitting in the optical/near-IR range. 
In contrast with the previous sample, this catalog includes also 
passive galaxies, and this explains the larger extention of the 
SFR-M, relation towards low star formation rates (blue trian- 
gles in fig. [2 1. These data lie on the predicted confidence region 
represented by the contour plot, however the bottom-left part of 
the SFR-M, diagram remains unsampled due to observational 
incompleteness. This observational limitation does not allow 
us to test one of the striking feature of the SFR-M, relation 
predicted by hierarchical models of galaxy formation, namely, 
the behaviour of the scatter. In fact, its increase with decreas- 
ing stellar mass is due to the large variety of star formation 
histories corresponding to the low-mass population, while the 
similarity of the star formation histories of high-mass galaxies 
results in smaller scatter in the SFR-M* plane. This is espe- 
cially true at high redshifts z > 6; at such epoch, the small 
number of progenitors of present-day low-mass galaxies (and 
their different collapse time) yield a large variance in the star 
formation histories, while the larger number of progenitors of 
high-mass galaxies (already collapsed since they formed in bi- 
ased regions of the density field) result in a smaller variance 
among the different star formation histories (see fig. [TJ. 
The dashed and long-dashed lines in figure |2| indicate the best- 
fit of the galaxy main sequence obtained by Rodigh iero et al.| 
( |2011| l and |Santini et al.| ( |200 9) respectively, while the solid 
line shows the spline fit to the peaks of the SFR (SSFR) dis- 
tributions of model galaxies as a function of stellar mass (log 
SFR=0.8*logM,-7.5). The slopes of the observed and predicted 
SFR-M, relations are similar (0.8 < a <0.9), while the nor- 



malization of the predicted relation is a factor of 4-5 lower 
than the observed relations. The offset between the predicted 
and observed SFR-M, relations in the local universe (z - 0.1 
Brinchmann et al.|2004[ l reduces by a factor of about two, im- 



plying that the evolution of the SSFR with redshift predicted 
by the model is slower than that observed at z < 2. The local 
distribution of model galaxies in the SFR-M, plane shows a bi- 
modal behaviour, with a clear separation between actively star 
forming and passively evolving galaxies, as seen in observa- 
tions (fBrinchma nn et al.T20 04). 

Similar offsets between the observed SFR-M, relation and 
those expected by various kind of SAMs and hydrodinamical 
simulations has been previously reported in the literature (see 
e.g. |Daddi et al.|2007||Dave|2008||Fontanot et al. 12009} |Damen 
et al.|2009||Santini et al.|2009[|Liii et al.|2012| >. Understanding 



whether this mismatch originates from systematic effects in the 
stellar mass and/or SFR determinations, or from incomplete- 
ness in our basic picture of galaxy assembly is a very difficult 
task. The amplitude of this offset depends on the the technique 
used to estimate the SFR and stellar mass and on the sample 
selection (see the offset between the dashed and long-dashed 
lines). In our analysis of the role of starbursts relative to qui- 
escently starforming galaxies (section |3~2| we shall normalize 
both the model and observed SFRs to their main sequence val- 
ues, and we will address this mismatch in more detail in section 

El 

To understand the role of starbursts and quiescently starform- 
ing galaxies in deriving the SFR-M, relation we separately 
show in figure [3] the relations obtained selecting from the 
Monte Carlo simulations galaxies dominated by the quiescent 
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Fig. 3. The predicted SFR-M* relation at 1.5 < z < 2.5 obtained selecting from the model galaxies dominated by the quiescent 
(left) and burst (right) mode of star formation. Contours and solid line as in figure [2] 



mode of star formation (SFR q > SFRb, left panel), and galaxies 
dominated by the burst mode of star formation (SFRb > SFR q , 
right panel). The former show a tight correlation between the 
SFR and the stellar mass. This is a natural outcome, being this 
component a fairly steady function of time (see the smooth 
component of the SFHs in fig. [TJ. The scatter of this relation 
originates from the distribution of cold gas at a given stellar 
mass, and ultimately stems from the stochastic nature of the 
merging histories of galaxies. When we select galaxies where 
SFRb > SFR q we obtain a more scattered distribution in the 
SFR-M* plane. In these galaxies the star formation rate is not 
only determined by the amount of cold gas available in the 
galaxy disk, but it is also regulated by the amount of such 
gas which is destabilized during galaxy interactions (see eq.|2|. 
Minor mergers and fly-by events, which dominate the statistics 
of encounters in the hierarchical clustering scenario, induce a 
lower fraction of destabilized gas compared to major mergers 
(note the dependence of / on the mass ratio m' /m in eq. [3J. 
This analysis show that the main sequence is mainly populated 
by quiescently starforming galaxies, while the loci well above 
and below it by galaxies experiencing major mergers and minor 
mergers/fly-by events respectively. 



3.2. The contribution of starbursts and quiescently 
starforming galaxies to the cosmic SFRD 

We now quantify the results shown in the previous section, 
and derive the relative contribution of starbursts and quies- 
cently starforming galaxies to the cosmic SFRD and comov- 
ing number density in the redshift interval 1.5 < z < 2.5. 
To perform a quantitative comparison we first need to define 



the sample of model starforming galaxies. Since the observa- 
tional results we compare with are derived from a combina- 
tion of samples adopting different criteria to select starforming 
galaxies (see Sect, 



3.1 



we have performed two different se- 
lections of model galaxies in order to reproduce, at least par- 
tially, the observational selection effects. The first corresponds 
to model galaxies with SSFR> SSFRms-2<x (see dashed line in 
fig. [2J3), where SSFRms is the value obtained for model galax- 
ies and <j = 0.3 dex is the scatter of the observed SFR-M* 
relation ( |Rodighiero et al.pOTTj |Daddi et al.|2004"| >. The sec- 
ond corresponds to model galaxies with SSFR> f#(z) (see 
the thin solid line in fig. ^> which indicates the correspond- 
ing SSFR value at z = 2), which select galaxies whose current 
SFR is stronger than its average past SFR (see e.g. |Fontana| 



et al. |2009] l . Within such samples, we derive the distributions 



of SSFR normalized to the value of main sequence galax- 
ies (SSFRms), in four equally spaced stellar mass bins from 
M„ = 10 1() M o to M, = 10 1L5 M Q (see fig. According to 



Rodighiero et al. (201 1 1, the galaxies observationally classified 



as starburst are those with SSFR^g > 4x SSFRms (vertical dot- 
ted lines in fig. [4). We then proceed to assess how galaxies ob- 
servationally classified as starburst compare with model galax- 
ies. To this aim, we derive separately the distributions of model 
galaxies dominated by the quiescent mode of star formation 
(blue histograms) and of galaxies dominated by the burst mode 
of star formation (red histograms). We find that the threshold 
SSFR SB > 4x SSFRms used by |Rodighiero et al.|p0TT) is in- 



deed effective in filtering out quiescent galaxies, independently 
of the criterion adopted to define our sample of starforming 
galaxies; in fact, in the model SSFR values above this threshold 
are obtained only through the burst component of the star for- 
mation. However this criterion fails to select a considerable part 
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Fig. 4. SSFR distributions of model galaxies, at 1.5 < z < 2.5, 
in four mass bins. The blue histograms refer to galaxies dom- 
inated by the quiescent mode of star formation, the red his- 
tograms to galaxies dominated by the burst mode of star forma- 
tion, and the green histograms to galaxies with an interaction- 
driven star formation component. The solid vertical lines show 
the position of the main sequence, while the vertical dotted 
lines show the SSFR threshold used to identify starburst galax- 
ies: SSFRsb > 4x SSFR M s. The solid distributions are obtained 
selecting galaxies with SSFR> ?h(z) _1 , while dashed distribu- 
tions are obtained selecting galaxies with SSFR> SSFRms-0.6 
dex. 



of the burst-dominated starforming galaxies where the galaxy 
interactions do not boost the SFRs above the quiescent values 
(minor mergers and fly-by events). The above results enlighten 
the physical difference between galaxies observationally clas- 
sified as starbursts and galaxies with an interaction-driven star 
formation component. The latter also include galaxies not dom- 
inated by the burst mode of star formation, as it is indicated 
by the distributions of galaxies with a burst component of star 
formation (SFRb > 0, green histograms in fig. [4|. In the fol- 
lowing we will derive the relative contribution to the comoving 
number density and to the cosmic SFRD in the redshift range 
1.5 < z < 2.5 for all these classes of starforming galaxies: (i) 
main sequence galaxies (SFR q > SFRb); (ii) starburst galax- 
ies (SSFRsb > 4x SSFR M s), (iii) burst-dominated starforming 
galaxies (SFRb > SFR q ); and (iv) galaxies with an interaction- 
driven star formation component (SFRb > 0). 
The predicted relative contribution of the above classes of 
galaxies to the comoving number density of starforming galax- 
ies is shown in figure [5^. The fraction of starbursts (shaded 
region) and burst-dominated galaxies (hatched region) with re- 
spect to the total number of starforming galaxies is compared 
with the observational results by Rodighiero et al. ( 2011) . 
While the predicted fraction of starburst galaxies remains 
around Nsb/Nms+sb ~ 1%, the exact values predicted by the 
model depend on the criteria adopted to define the global sam- 



ple of starforming galaxies (i.e., N M s+sbX defining the upper 
and the lower envelope of both the shaded and the hatched re- 
gions. The smaller fraction of starburst galaxies obtained with 
the selection SSFR> SSFRms-2<x compared to the selection 
SSFR> tnizT 1 is due to the larger number of "on-sequence" 
galaxies selected in the first case (especially for large masses, 
see solid and dashed histograms in fig. [4}, while the number of 
galaxies N5 b remains almost unchanged. 
The figure shows that, even in a hierarchical model connecting 
the physics of galaxies to their merging histories, the fraction 
of galaxies classified as starburst remains within values ~ 1%, 
even lower than the observational values (~2%); values around 
2-3% are obtained when burst-dominated galaxies are consid- 
ered. This is due to a twofold reason: i) in hierarchical models, 
the statistics of merging events is dominated by minor episodes 
which do not produce SSFR above the threshold adopted in the 
observations; ii) the duration of bursts t as 10 7 - 10 8 yrs is 
much smaller than the typical time scale between two encoun- 
ters, so that the probability to find a galaxy in a starburst phase 
is correspondingly small. As for former point i), we show in fig- 
ure [5^ the fraction of the galaxy population with a star forma- 
tion component induced by galaxy interactions, SFRb > (tri- 
angles). When all galaxy interactions are considered, includ- 
ing minor merging and fly-by, the fraction of galaxies with an 
interaction-driven star formation component increases with M* 
to reach values w 10 _1 . The trend with M, is due to the larger 
cross section for interactions of massive galaxies (see eq. [4]) 
and to the fact that massive galaxies are found in dense envi- 
ronments with enhanced galaxy interaction rate. The merging 
histories of massive galaxies are dominated by minor merging 
events which are less effective in destabilizing the gas in the 
galaxy disks (see eq. [3]), therefore a lower number of massive 
galaxies satisfy the SFRb > SFR q criterion, and this explains 
the flattening of the hatched and shaded regions with increasing 
stellar mass in figure [5^. 

As for the starburst duty cycle, the above value ~ 1CT 1 
of the fraction of galaxies with SFRb > implies a value 
Tb/T r ~ 0.1 for the ratio of the burst duration T& over the 
interaction time scale r r (see sect. 2.2). This is indeed what is 
expected from both theoretical arguments and observational 
hints. In fact, the interaction time scale t> = (n'LVy 1 can 
be estimated assuming a simple geometrical cross section 
X * nr 2 . For typical values of the galaxy gravitational radius 
r ~ 0.1 Mpc and of the galaxy relative velocities V ~ 100 
km/s, and assuming a density of a few galaxies per Mpc 3 , 
we obtain an estimate r r ~ 10 9 yr. When compared with the 
galaxy crossing time (determining the duration of bursts, see 
sect. 2.2) T[, ~ 10 8 yr, a value t^/tv m 0.1 is obtained for the 
duty cycle, as expected. This value is also consistent with the 
duty cycle of A GN measured at z ~ 1-2 with lumin osity 



L > 10 4J erg/s (Martini et al. 



2009 



Eastman et al. 



2007 1, as 



expected in our framework where both starbursts and AGNs 
are triggered by galaxy interactions. 

We now proceed to estimate the contribution of the different 
classes of starforming galaxies to the cosmic SFRD at 1.5 < 
z < 2.5, shown in figure [5]?. The interaction-driven scenario 
predicts that the contribution of the entire interacting galaxy 
population is ~15% irrespective of stellar mass. At lower stellar 
masses this quantity is mainly contributed by burst-dominated 
galaxies, while at higher stellar masses the contribution of the 
quiescent mode of star formation becomes important as it is 
indicated by the decrease of the ratio SFRDsb/SFRDms+sb 
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Fig. 5. Left: fraction of starburst (shaded region), burst-dominated (hatched region), and interacting (triangles) galaxies as func- 
tion of the stellar mass at 1.5 < z < 2.5. Right: contribution to the cosmic SFRD of the same classes of starforming galaxies. 
The lines indicate the model predictions adopting the SSFR> tnizY 1 (solid lines) and the SSFR > SSFRms-2o" (dashed lines) 
criterion to define the global sample of starforming galaxies. The data points indicate the results obtained by Rodighier o et al.| 
( 20TT) for starburst galaxies (SSFR>4x SSFR MS ). 



with increasing M, for the SFR b > SFR q and SSFR SB > 4x 
SSFRms populations. In fact, for the burst-dominated starform- 
ing galaxy population the contribution to the cosmic SFRD is 
-15% for M, ~ IO^Mq and (6-10)% for M» ~ 1O U ' 5 M , 
and for starburst galaxies it is slightly lower: ~ 12% for M, ~ 
1O 1O M and (4 -6)% for M. ~ 10 1L5 M Q . Therefore, this sce- 
nario is consistent with the observational finding that the cos- 
mic SFRD in the redshift range 1.5< z <2.5 is mainly con- 



tributed by quiescently starforming galaxies (Rodighiero et al. 
201 1) . However, it must to be noted that, given the lower frac 



tion of starburst galaxies predicted by the model, the match 
between the predicted and observed starburst contributions is 
determined by the lower value of the SFR of model main se- 
quence galaxies (see figure[2]). 

Since the burst mode of star formation is intimately related to 
the galaxy interaction rate, which is expected to increase with 
redshift, we investigate how these contributions change as a 
function of the cosmic epoch. To address this issue, we derive 
separately the contribution to the cosmic SFRD of quiescent- 
dominated and burst-dominated systems as a function of red- 
shift (dashed and dotted lines in figure ^respectively). We find 
that quiescent systems dominate the global SFRD at all red- 
shifts. However, the contribution of the burst-dominated popu- 
lation increases with redshift, rising from < 5% at low redshift 
(z < 0.1) to ~ 20% at z Z 5. This is a typical feature of the 
hierarchical clustering scenario where the starburst events are 
triggered by galaxy interaction during their merging histories; a 
similar evolution was obtained by Hopkins et al. (2010 1 on the 
basis of cosmological hydrodinamical simulations. The phys- 



ical origin of this behaviour can be understood as follows: at 
high redshift both the interaction rate r" 1 (eq.|4]) and the frac- 
tion of destabilized gas / (eq. [3]) are large; the former due to 
the large densities, the latter due to the large ratio m' /m ~ 1 
characteristic of this early phase (z > 3) when galaxy inter- 
actions mainly involve partners with comparable mass (major 
mergers). At lower z the decline of the interaction rate and of 
the destabilized gas fraction suppresses the burst mode of star 
formation thus lowering its contribution to the total SFRD. 
An immediate implication of the above is that the fraction 
of the stellar mass formed through the SFR associated with 
the destabilized cold gas during galaxy interactions increases 
with redshift. This is illustrated in figure [7] where the con- 
tours show the average values of the ratio Mb urs t(z)/M»(z) as 
a function of M, and z- The mass Mb ulst (z) is the stellar mass 
that a galaxy has formed through the burst mode of star for- 
mation up to that redshift, and M,(z) is the total galaxy stel- 
lar mass at that z. From this figure one can infer that ~10% 
of the final (z = 0) galaxy stellar mass has been formed dur- 
ing bursts for 10 10 M o <M, < 10 115 M G . At higher redshift 
(z > 4) this fraction increases up to ~20%. To illustrate how 
the ratio Mb U rst(z)/M t (z) varies during the formation of galax- 
ies with different final masses, we also show the growth histo- 
ries of typical galaxies with final mass: M*(z = 0) = 10 12 M o , 
M„(z = 0) = 10 n M G , and M t (z = 0) = 10 ll) M Q (red, blue, and 
purple lines in fig. [7}. These paths illustrate that at high redshift 
(z > 2) and for massive objects large values of Mb ur st(z)/M»(z) 
are expected, since such massive galaxies formed in biased re- 
gions of the density field where high-redshift interactions are 
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Fig. 6. Contribution to the cosmic SFRD of galaxies with 
SFR>1.4xlO" 6 M G /yr (corresponding to UV l uminosity L rrv > 
10 15 W Hz 1 adopting the calibration given in Kennicutt 1998 1 
as a function of redshift (solid line). Dashed line indicates the 
contribution of galaxies dominated by the quiescent mode of 
star formation, while dotted line indicates the contribution of 
galaxies dominated by the burst mode of star formation. 



extremely effective in triggering starbursts. For less massive 
galaxies the high-redshift value of the Mb ulst (z)/M H ,(z) ratio 
progressively lowers, since these galaxies formed in less dense 
environments. 



4. Discussion 

In this section we discuss the origin of the mismatch between 
the normalization of the predicted and observed galaxy main 
sequences which was highlighted in section 3.1 Similar mis- 
match between the observed SFRs and those expected by var- 
ious kind of SAMs and hydrodinamical simulations has been 

< 2 (see e.g. |Daddi et aL 



also reported by several authors at z 

200 7] |Dave|l2008l IFontanot et al |[2009] |Damen et al.||20091 
Santini et al'|2009| Lin et al.|2012| i. This discrepancy either in- 
dicates an incompleteness in our knowledge of galaxy assem- 
bly, or systematic effects in stellar mass and/or SFR determi- 
nations, or both. 

On the observational side, the uncertainties in stellar mass and 
SFR estimates may be responsible for one of the main puzzles 
that appear in present-day observational cosmology: the mis- 
match between the observed stellar mass density and the inte- 
grated SFRD (e.g.|Hopkins & Bea com 2006] |Fardal et al.|2007| 
Wilkins et al. 2008 ; San tini et al.|2Q12p . In principle, these two 
observables represent independent approaches to studying the 
mass assembly history from different points of view. However, 
the integrated SFRD, after considering the gas recycled frac- 
tion into the interstellar medium, appears to be higher than the 
observed stellar mass density by a factor of about 2-3 at z < 2 
(San tini et al.|2012] >. Moreover, if the merging contribution to 
the stellar mass build-up is accounted for (Drory & Alvarez 
2008) the agreement gets even worse. Intriguingly, the inte- 
grated SFRD exceeds the observed stellar mass density, imply- 
ing that either the SFR are overestimated, or the stellar masses 
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Fig. 7. The predicted average values of the ratio 
M[, urs t(z)/M t (z) as a function of redshift and galaxy stel- 
lar mass. The four filled contours correspond to equally 
spaced values of the ratio Mb ulst (z)/M,(z): from 0.05 for the 
lightest filled region to 0.2 for the darkest. The lines illustrate 
the growth histories of a typical galaxy with final mass: 
M»(z = 0) = 10 12 M o (red), M»(z = 0) = 10 n M o (blue), and 
M,(z = 0) = 10 10 M G (purple). 



are underestimated. Both effects would reduce the offset be- 
tween the observed and predicted SFR-M* relations at least by 
a factor of 2. 

On the theoretical side, we have investigated whether the mis- 
match can be ascribed to the implementation of the star forma- 
tion law usually adopted in SAMs. In this respect, the most crit- 
ical point is the assumption that the star formation law in eq. (fT|) 
- derived from the local observed Kennicutt relation (Ken nicutt| 
1998) - is valid at any redshift. To investigate the effect of re- 
laxing this assumption, we have varied the high-redshift star 
formation timescale (the free parameter q in eq.[T| so as to ob- 
tain star formation efficiencies SFR/m c above the local value 
for z > 2. However, even increasing the efficiency up to a fac- 
tor of 4, we do not find an appreciable up-shift of the model 
main sequence. In fact, increasing the high-z star formation re- 
sults into an earlier consumption of the the cold gas reservoir 
of model galaxies, an effect which balances the increased star 
formation efficiency leaving the normalization of the SFR-M* 
relation almost unchanged. 

It is also interesting to examine the impact of other star forma- 
tion laws on the shape and normalization of the main sequence. 



This was investigated by |Lagos et al. ( 201 lb} by implementing 
the empirical relatio n of |Blitz & Rosolowsky| ( |2006[ ) and the 
theoretical model o f rKrumholz et al.| l|20097 in the SAMs of 
Baugh et al. (2005 1 and |Bower et al.| (|2006). These star forma- 
tion laws are of the form Ssfr = vsF^moi = vsFfmol^gas> where 
Ssfr and E gas are the surface densities of the SFR and total cold 
gas mass, vsf is the inverse of the star formation time scale 
for the molecular gas, and f mo i = £ mo //£ gfls is the molecular- 
to-total gas mass surface density ratio which depends on the 
hydrostatic gas pressure ( |Blitz & Rosolows ky 2006 ) or on the 
balance between the dissociation of molecules due to the inter- 
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stellar far-UV radiation and their formation on the surface of 
dust grain ( |Krum holz et al. 2009 ). They found that the main se- 
quence appears insensitive to the star formation law. Thus, tun- 
ing the quiescent star formation law does not seem to provide a 
solution to the mismatch between the normalization of the pre- 
dicted and observed galaxy main sequence. It is however possi- 
ble that other ways of modelling starburts, e. g. starburst driven 
by disk instabilities (e. g.|Bower et al.|2 006 ; Lag os et al.|2011b| 



Guo et al.|20"TT]|Fanidakis et al.|2012[|201 lHHirschmann et al] 
2012), could have an impact on the shape and normalization 
of the main sequence (see [Lagos et a l.|201 lb) . While the lat- 
ter constitutes an interesting complement to the merger-driven 
starbursts, its consistence with the observed fraction of star- 
bursts is still to be investigated. 



It is also to be noted that the offset between the observed 
and predicted galaxy main sequence can be due not only to an 
underestimation of the predicted SFRs but also to an overesti- 
mation of the predicted stellar masses. A physical process that 
can reduce the stellar mass of galaxies is the disruption of stel- 
lar material from merging satellites due to tidal stripping (e.g. 
Henriques & Thomas 2010). Evidence for the importance of 
this process in galaxy formation comes from the existence of a 
diffuse population of intra-cluster stars ( Zwicky|1951 Durrell 
et al.| 2002] ). Indeed, gas-dynamical simulations (e.g. Moore 
et al.|| 1996 1 generally agree that these stars have been contin- 



ually stripped from member galaxies throughout the lifetime 
of a cluster, or have been ejected into intergalactic space by 
merging galaxy groups , rather than having formed in the intra- 
cluster medium. 



An alternative possibility is that the cold gas fractions 
predicted by the models at z ~ 2 underestimates the real 
values, due to a too early conversion of cold gas into stars 
at earlier cosmic times. However, limiting the cooling and 
star formation efficiency at high redshifts is an extremely 
challenging task in the hierarchical clustering scenarios. 
These predict the high redshifts progenitor galaxies to be 
characterized by low virial temperature and high densities, 
making the gas cooling extremely efficient at high redshifts 
z > 3; the short dynamical timescales and the frequent merging 
events rapidly convert such available cold gas into stars, 
leaving only a residual fraction of galactic baryons in the form 
of cooled gas available for star formation at z ~ 2. Although 
such an over-cooling problem (already pointed out in the 
seminal paper by White & Rees 1978 ) is alleviated by stellar 
feedback expelling part of the cold gas from the starforming 
galactic regions, the efficiency of the latter process is severely 
limited at high redshifts due to the compactness of dark matter 
haloes resulting in large escape velocities. In fact, such a 
long-standing problem of the hierarchical scenario of galaxy 
formation leads to the over-prediction of low-mass galaxies 
at z > 1.5, as shown by the comparison of our and other 
SAMs with the observed stellar mass function (Fontana et al. 



2006; Fontanot et al. 2009; Marchesini et al. 2009; Guo et al. 



201 l] |Santini et al.|2012| l. The excess of starforming low-mass 
galaxies at z > 1-5 reflects into an excess of red low -mass 
galaxies at z ~ 0, even when a large feedback is implemented 
in the models (Gu o et al.|2011) . 

Assessing whether the evolution of the cold gas in galaxies 
is the key to explain (at least partially) the SAMs under- 
prediction of SSFRs at z < 2 would require a comparison 
between the predicted and the observed galaxy cold gas 
fraction as a function of redshift. The evolution of the cold 
gas content of galaxies represents also one way of testing 
time-changing star formation laws, since the latter determine 
the rate at which gas is converted into stars. An increase in 
the star formation time scale leads to larger cold gas masses 
in disks and vice versa (see e.g. Lagos et al. 201 la|b 



Fu 



|et al.| |2012| i. In models where starbursts are triggered by 
disk instabilities, a longer star formation time scale in disks 
can lead to lower final gas content in some cases, being the 
disk more prone to instabilities if the total mass is larger 
(Lagos et al. 2011b). The predictions for such quantity in 
our model are shown in figure [8] in terms of the gas fraction 
parameter f gas = M gas /(M gas + M„), for model galaxies with 
SSFR> t H (zY l and M, > 10 I0 M Q , as predictions to be tested 
with future observations. 
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Fig. 8. Gas fraction distributions as a function of redshift of 
model galaxies with SSFR> t H {z)~ l and M, > 10 10 M o . The 
filled contours correspond to equally spaced values of the den- 
sity (per Mpc 3 ) of galaxies in logarithmic scale: from 10" 3 
for the lightest filled region to 10~ 21 for the darkest. The 
plotted value of log f gas is the median value for sources 
in four redshift bins: log f gas =-036+ f 4 (0< z <0.5); log 
/ gflS =-0.32!^ (0.5< z <1); log / gflS =-0.29^ (1< z <1.5); 

log/ gfls =-0. 25^05 (1.5< z <2.5). The vertical error bars show 
the interpercentile range containing 68% of the sources. 



5. Conclusions 

We have investigated the relative importance of the quiescent 
mode and burst mode of star formation in determining the SFR- 
M, relation and the evolution of the cosmic SFRD, predicted 
under the assumption that starburst events are triggered by 
galaxy encounters during their merging histories. The latter are 
described through Monte Carlo realizations, and are connected 
to gas processes and star formation using a SAM of galaxy for- 
mation in a cosmological framework. The main results of this 
paper are: 
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- Hierarchical clustering models reproduce the slope and the 
scatter of the SFR-M* relation, however, they under-predict 
the normalization by a factor of 4-5 a z — 2 (Daddi et al. 
20071 |Dave|2008HFontanot et al.|2 009 ; Da men et al.|2QQ9[ 



Santi ni et al.||2009| |Lin et al.||2012| i A possible theoreti- 



cal explanation of this mismatch is that the amount of cold 
gas in galaxy disks predicted by these models at z — 2 un- 
derestimate the real values. We derive the prediction for 
the evolution of the gas fraction distribution of starforming 
galaxies, as prediction to be tested with future observations. 

- The tight correlation between SFR and M, of the galaxies 
on the main sequence is determined by the quiescent com- 
ponent of the star formation. Galaxies dominated by the 
star formation component induced by galaxy interactions 
populate the regions well above (major mergers) and below 
(minor mergers and fly-by events) the sequence. 

- The predicted SSFR distributions of starforming galax- 
ies indicate that galaxies observationally classified as star- 
burst on the basis of their distance from the main se- 
quence (SSFRsb > 4 x SSFRmsX are indeed dominated by 
the burst component of the star formation. However, this 
criterion fails to select part of the burst-dominated star- 
forming galaxies where galaxy interactions do not strongly 
boost the star formation rate. 

- Hierarchical clustering scenarios, connecting the proper- 
ties of galaxies to their merging histories, naturally yield 
fraction of starburst galaxies at z ~ 2 of ~1%, in agree- 
ment with the observational estimate (2% Rodighiero et al. 
|201 1[ ). This low fraction is due to the fact that in this sce- 
nario, the statistic of merging events is dominated by minor 
episodes which do not produce SSFR above the observa- 
tional threshold, and to the small probability of finding a 
galaxy in a starburst phase due to the small value of the du- 
ration of the burst (f ^ 10 7 - 10 8 ) yrs compared with the 
typical time scale between galaxy encounters. 

- Quiescently starforming galaxies dominate the global 
SFRD at all redshifts. The contribution of the burst- 
dominated systems increases with redshift, rising from < 
5 % at z < 0.1 to ~ 20% at z > 5. This behaviour is deter- 
mined by the increase of the galaxy interaction rates and of 
the effectiveness of galaxy interaction in destabilizing the 
cold gas at high redshift. The fraction of the final (z = 0) 
galaxy stellar mass which is formed through the burst com- 
ponent of star formation is -10% for 10 10 M o <M„ < 
10 1L5 M G . 

These findings do not imply that galaxy mergers/interactions 
play a lesser role in the formation of stars in galaxies. Indeed, 
for each galaxy the star formation (both the quiescent and burst 
component) is driven by the cooling rate of the hot gas and by 
the rate of refueling of the cold gas, which in turn is intimately 
related to the galaxy merging histories. 

We found that the interaction-driven model for starburst galax- 
ies provides strict predictions for the evolution and the mass 
dependence of the starburst contribution to the cosmic SFRD, 
which represent effective tools to test this scenario with fu- 
ture observations of large and complete samples of starform- 
ing galaxies at z > 2. To reach this goal, it is also necessary 
a reliable observational diagnostic which is able to distinguish 
between galaxies dominated by the quiescent and burst mode 
of star formation. The measure of the star formation efficiency 
and of the IR star formation compactness could represent eligi- 
ble candidates for this purpose ( Daddi et al.|2010||Genzel et al. 
|2010 ; Elbaz et al.||20rT) . The first is based on measurements 



of the mass and density of molecular gas at high redshift, and 
on the poor knowledge of the CO luminosity to H2 conversion 
factor. The technique of separating starbursts and quiescently 
starforming galaxies based on their IR star formation compact- 
ness seems to be promising, since ALMA will provide power- 
ful tools to measure the spatial distribution of star formation in 
distant galaxies at high angular resolution, making it possible 
to estimate the compactness of the star formation regions. 
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